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(N ■ Abstract 

in 

^ ■ A study of b semileptonic decays into D, Dn^ and Z^*7r ± final states is presented. 

The D°, D + and D* + mesons are exclusively reconstructed in Z decay data 
recorded from 1992 to 1995 in the DELPHI experiment at LEP. The overall 
branching fractions are measured to be: 

X ■ 

V ; BR(b -> D £~u e X) = (7.04 ± 0.34 (stat) ± 0.36 (syst.exp) ± 0.17 (BR D ))% 

BR(6 -> D + £~i? i X) = (2.72 ± 0.19 (stat) ± 0.16 (syst.exp) ± 0.18 (BR D ))% 

43 ' 



BR(6 D* + n? t X) = (2.75 ± 0.17 (stat) ± 0.13 (syst.exp) ± 0.09 (BR D ))% 



where the .D and D + results include also contributions from D*° and D*^ 



decays. A fit to the distribution of the 7r ± impact parameter to the primary 
interaction vertex provides a measurement of the b semileptonic branching frac- 
tions into the D°tt ± X, D + ti ± X and D* + 7r ± X final states. Assuming that single 
pion decay modes of B mesons dominate, the partial rates for B — > D-ni"vi and 
B — > D*tt£~i7£ have been obtained, corresponding to a total branching fraction: 

BR(B -> Ditrvi) + BR(5 D*ir£-j7 e ) = (3.40 ± 0.52 (stat) ± 0.32 (syst))% . 

This result agrees well with the observed difference between the total B semilep- 
tonic branching fraction and the sum of the B — > Di~vi and D*i~vt branching 
fractions. 



(Phys. Lett. B475(2000)407) 
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1 Introduction 

The study of B meson semileptonic decays into any Dtt or D*n final state is interest- 
ing for several reasons. Present measurements of B semileptonic decays into D£~vi and 
D*£~i?e imply that these final states account for only 60% to 70% of all B semileptonic 
decays [|TJ . The remaining contribution could be attributed to the production of higher 
excited states or non-resonant D^tt final states, hereafter denoted D**. However the 
ALEPH measurement of the D**£~v^ branching fraction does not fully account for the 
observed discrepancy ||. The ratio of branching fractions of B decays into D* + 7T~£~i?e 
over all D* + £~X final states [] is also a significant contribution to the systematic uncer- 
tainty on r B o, Am d or V ch measurements [|l]J|. 

This paper describes a measurement of the branching fraction of B —> D**£~vi decays 
in the DELPHI experiment at LEP. The decays of the D°, D + and D* + into -D°7r+ are 
exclusively reconstructed 0. The analysis of D** — > D^ir^ in B semileptonic decays | 
relies on the impact parameter of the 7r** candidate, defined as its distance of closest 
approach to the reconstructed primary interaction vertex. A similar technique has been 
applied previously in ALEPH [||,^| and DELPHI ||. The single pion final states D°n + , 
D + 7i~ or D* + 7i~ , denoted "right" sign, are expected to dominate the decay widths. But 
pion pair emission, such as Dtc + 7t~ , is also allowed and could provide "wrong" sign D°tt~ , 
D + 7i + or D* + 7i + combinations. Similarly D s orbitally excited states can decay into D°K + 
or D*°K + which can be distinguished from D ( -* S)0 n + if the kaon is identified. 

The overall semileptonic branching fractions of a b quark into D°, D + or D* + final 
states are also presented in this paper. 



2 The DELPHI detector 

The DELPHI detector has been described in detail elsewhere |||7]; only the detectors 
relevant to the present analysis are briefly described in the following. The tracking of 
charged particles is accomplished in the barrel region with a set of cylindrical tracking 
detectors whose axis is oriented along the 1.23 T magnetic field and the direction of the 
beam. 

The Vertex Detector (VD) surrounds a Beryllium beam pipe with a radius of 5.5 cm. 
It consists of three concentric layers of silicon microstrip detectors at radii of 6, 9 and 
11 cm from the beam line. In 1991-1993 all the VD layers were single-sided with strips 
parallel to the beam direction. In 1994 and 1995, the innermost and the outermost layers 
were replaced by double-sided silicon microstrip modules, providing a single hit precision 
of about 8 /im in r<p, similar to that obtained previously, and between 10 /mi and 20 /xm 
in z H 0. For polar angles between 44° and 136°, a track crosses all the three VD layers. 
The innermost layer covers the polar angle region between 25° and 155°. For charged 
particle tracks with hits in all three rip VD layers, the impact parameter precision is ||: 

°r<j> = a 3/2a ® b C 1 ) 
p sin ' 

where a — 61 ± 1 /im, 6 = 20 ± 1 /mi and p is the momentum in GeV/c. 

throughout the paper charge-conjugate states are implicitly included; I indicates an electron or a muon, not the sum 
over these two leptons. 

2 7r* denotes the charged pion from the _D*+ — > D°7r+ decay. 

3 D stands for D° or -D+; 7r„* denotes the charged pion from the decay of a higher excited state of charmed meson or 
from a non-resonant D^n final state. 

4 In the DELPHI coordinate system: z is along the beam line, <f> is the azimuthal angle in the xy plane, transverse to 
the beam axis, r is the radius and 8 is the polar angle with respect to the z axis. 
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The Inner Detector is placed outside the VD between radii of 12 cm and 28 cm. It 
consists of a jet chamber giving up to 24 spatial measurements and a trigger chamber 
providing a measurement of the z coordinate. The VD and ID are surrounded by the 
main DELPHI tracking chamber, the Time Projection Chamber (TPC), which provides 
up to 16 space points between radii of 30 cm and 122 cm. The Outer Detector (OD) at 
a radius of 198 cm to 206 cm consists of five layers of drift cells. The average momentum 
resolution of the tracking system is a(l/p) < 1.5 x 10~ 3 (GeV/c) -1 in the polar angle 
region between 25° and 155°. The tracking in the forward (11° < 9 < 33°) and backward 
(147° < 9 < 169°) regions is improved by two pairs of Forward drift Chambers (FCA and 
FCB) in the end-caps. 

Hadrons were identified using the specific ionization (dE/dx) in the TPC and the 
Cherenkov radiation in the barrel Ring Imaging CHerenkov detector (RICH) placed be- 
tween the TPC and the OD detectors. 

The muon identification relied mainly on the muon chambers, a set of drift cham- 
bers giving three-dimensional information situated at the periphery of DELPHI after 
approximately 1 m of iron. 

Electron identification relied mainly on the electromagnetic calorimeter in the barrel 
region (High density Projection Chamber HPC) which is a sampling device having a 
relative energy resolution of ±5.5% for electrons with 45.6 GeV/c momentum, and a 
spatial resolution along the beam axis of ±2 mm. 



3 Event selection and simulation 

Charged particles were required to have a measured momentum between 0.3 GeV/c 
and 50 GeV/c, a relative error on momentum less than 100%, a track length in the TPC 
larger than 30 cm and a distance of closest approach to the interaction point of less than 

4 cm in r and less than 10 cm in z. 

Hadronic events were required to have at least five charged particles with momentum 
greater than 0.4 GeV/c and a total energy of the charged particles (assumed to be pions) 
greater than 12% of the collision energy. A total of Nz = 3.51 million hadronic events 
was obtained from the 1992-1995 data. Simulated hadronic events were generated using 
the JETSET 7.3 Parton Shower program 8.5 million Z — > qq and 4.0 million Z — > 



bb generated events, corresponding to seven times the available statistics in real data for 
bb final states. The B meson mean lifetime was set to t^ c = 1.6 ps. The generated events 
were followed through a detailed detector simulation and then processed through the 
same analysis chain as the real data. The hadronic event selection efficiency was thus 
estimated to be tz = 95.7%. The data sample contained also 0.2% of r pair events and 
0.2% of Bhabha events. 

The primary interaction vertex was computed in space for each event using an iterative 
procedure based on the x 2 of the fit. The average transverse position of the interaction 
point, known for each fill, was included as a constraint during the primary vertex fit. 
In order to increase the bb purity of the selected sample, using the impact parameter of 
all measured charged particle tracks in the event, the probability that all these tracks 



originate from the primary vertex was required to be smaller than 0.1 [IT|. This selection 
retains 15% of Z — > uu, dd and ss events, 48% of Z — > cc events and 94% of Z — > bb 
events. 

In order to estimate the reconstruction efficiencies and the invariant mass resolu- 
tions, dedicated samples of events containing a B meson decaying into D tt + £~uiX , 
D + it~ £~VfX or D* + n~ t~i>gX were generated. Physical backgrounds have also to be stud- 
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ied. These can be due to b — > cW~ decays followed by W~ — > cs and c — > £~z/^X (hereafter 
denoted 6 — > c — > £ background), or followed by c — ► £ + z/^X with — > DX (denoted 
6 — > c — ► 2 background). For this purpose, some dedicated samples of B — > DD S X or 
DDKX decays, with the D s or D decaying semileptonically, were generated. 

4 D(*k~ selection 

4.1 Lepton selection and identification 

Both muon and electron candidates were selected with a momentum larger than 
2 GeV/c. The lepton candidate was required to have at least one hit associated in 
the Vertex Detector. 

The muon identification algorithm is described in reference |7]] . A "loose" selection cri- 
terion provided an identification efficiency of (90 ±2)% for a probability of misidentifying 
a charged hadron as a muon of 1.2% within the acceptance of the muon chambers. 

A neural network procedure, combining information from several detectors, has been 
developed for electron identification. Electrons were identified with an efficiency of (65 ± 
2)% and a misidentification probability that a hadron be identified as electron of about 
0.4% @. 

The lepton transverse momentum relative to the meson momentum vector (as de- 
fined below) was required to be larger than 0.7 GeV/c. This cut reduced the contamina- 
tion of leptons from b semileptonic decay into r and from b^c^ioib^c^l decays. 

4.2 £)W decay channels 

The meson candidates were reconstructed in the following decay channels: D° — > 
K~7i + or K~7t + 7t + 7t~ (for D° not coming from a D* + decay), D + — > K~7t + tt + and 
D* + — ► D°7T+ with a D° decaying into K~n + , K~7t + 7t + tt~ or K~7t + (tt°) where the ir° 
was not reconstructed. In order to optimize the statistical precision of the measured 
production rates, slightly different selection criteria, as described below, were chosen in 
each £)W meson sample. 

Only charged particles with momentum vectors in the hemisphere defined by the lepton 
direction were considered for the reconstruction of charmed mesons. The kaon candidate 
from the D decay was required to have the same charge as the identified lepton. The kaon 
momentum was required to be larger than 1 (2) GeV/c in the D° (D + ) channel. The 
momentum of each pion from the D°/ + decay had to be larger than 1 GeV/c, except for 
the K~7r + 7r + 7r~ final state where the minimum momentum of candidate pions was lowered 
to 0.5 (0.3) GeV/c in the D° (D* + ) analysis. Any charged particle with a momentum 
between 0.3 GeV/c and 4.5 GeV/c and a charge opposite to that of the kaon was used 
as pion candidate for the D* + — > D°7r+ decay channel. 

To reduce the combinatorial background for all channels, except in the D* + —> 
(X _ 7T + )7r+ decay, the kaon candidate of the D was required to be identified according to 
the RICH and dE/dx information 0. In the D° -> K~7r + and D° -> J\-7r + (7r°) decay 
channels, the angle 9* between the K~ir + momentum vector and the kaon direction in 
the K~n + rest frame was required to satisfy the condition cos 6* > —0.8. For genuine 
D° candidates an isotropic distribution in cos 8* is expected whereas the background is 
strongly peaked in the backward direction. 

The 7T* candidate and at least two particles from the D 0//+ decay were required to have 
at least one hit associated in the Vertex Detector. 
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4.3 D vertex 

After the previous selections, a K~7t + , K~7t + tt + or K~7t + 7t + tt~ vertex was fitted in 
space. In the K~tt + tt + tt~ decay channel, in order to reduce the large combinatorial back- 
ground, the impact parameters of charged particle trajectories, relative to the common 
D° vertex, were required to be smaller than 150 /zm. In the D + — > K~7t + tt + channel, 
either these impact parameters had to be less than 100 /im or the \ 2 probability of the 
K~~7t + tt + vertex had to be larger than 10 -4 . 

The momentum vector of each particle, attached to the D vertex, was recomputed at 
this vertex. In each channel, the scaled D energy, Xe(D) = E(D)/Eb eam , was required 
to be larger than 0.15. 

The apparent decay length of the D° or D + candidate, AL, was computed in the 
plane transverse to the beam axis. It was given the same sign as the scalar product of 
the D momentum direction with the vector joining the primary to the D vertices. In the 
D* + channel, AL was required to be positive. In the D° and D + channels, which have 
a higher combinatorial background, the value of AL divided by its error was required to 
be larger than 1. 

4.4 B vertex 

Finally a D°£, D + £ or D°7r*£ vertex (denoted "B" vertex in the following) was fitted 
in space. The B decay length was defined, as above, as the signed distance between the 
primary vertex and the secondary D(n t )£ vertex in the plane transverse to the beam axis. 
This B decay length divided by its error was required to be larger than 1 for all channels. 
In order to reduce further the combinatorial background, the decay length divided by its 
error between the D and the B vertices was also computed: it was required to be larger 
than -1 in the D° samples and in the D* + — > (K~ n + -k + 71^)11+ sample, and to be positive 
in the D + sample. 

4.5 D invariant mass 

The selection of D* + £~ events relied on the small mass difference (AM) between 
the D* + and the candidate D°. On the contrary, D* + candidates were rejected from 
the D° and D + samples as follows: the D° candidates were rejected if at least one 
7T* particle was found in the event giving a AM value less than 160 MeV/c 2 ; in the 
D + — > K~7t + tt + sample, both K~ir + pairs were associated to the remaining ir + and a 
AM mass difference was computed, the K~tt + 7t + combination was rejected if at least 
one AM value was found smaller than 160 MeV/c 2 . 

Figure [I] shows the invariant mass (or mass difference) distributions in each of 
the previously selected D meson channels. In the D* + — > (L^7r + )7r+ channel, the 
K~7i + invariant mass had to be within 75 MeV/c 2 of the nominal D° mass. In the 
D* + — > (Lf~7r + (7r ))7r+ channel, the K~n + invariant mass had to be between 1500 and 
1700 MeV/c 2 . In the D* + -> (Af"7r + 7r + 7r")7r+ channel, the AM mass difference had 
to be within 2 MeV/c 2 of the nominal D* + — D° mass difference. The invariant mass 
of the D° — > K~tt + channel has a resolution of about 25 MeV/c 2 whereas it is less 
than 15 MeV/c 2 in the K~ ~7r + 7r + 7r final state. Thus, in the particular case of the 
D* + — > (K~tt + )tt+ decay channel, the K~tt + invariant mass was constrained to the 
D° mass value and a constrained (D°7r*£) kinematic fit was performed. This improved 
the resolution of the (D* + — D°) mass difference by 30% in this channel. 
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A clear signal corresponding to D£~ events is observed in each distribution (data 
points), whereas the wrong sign Di + combinations (hatched histograms) present a much 
smaller D meson contribution. The right sign invariant mass distributions were fitted 
with a signal component described by the sum of Gaussian functions, and a combinatorial 
background parameterised with a polynomial form. In the D° — ► Kir and Kirmr samples, 
the contribution from missing 7r° appears as a "satellite" peak for mass values smaller 
than the nominal D° mass. This contribution was parameterised as the sum of Gaussian 
functions with their parameters fixed according to the simulation. In each channel, the 
relative amounts and relative widths of the Gaussian functions describing the D signal 
were tuned according to the simulation. The free parameters of the fits were thus the 
coefficients of the polynomial background, the normalisation of the "satellite" peak (in 
the Kti and Kittttt invariant mass distributions), the average width and mean value of the 
signal shape and the number of D signal candidates. For each decay channel, the mass 
distributions of the wrong sign D£ + events were fitted with the same shape parameters 
as the right sign signals. This allowed the contribution of the fake lepton events to be 
determined and then subtracted. The observed numbers of D mesons, within the quoted 
range around the D°/ + mass and D* + — D° mass difference, is indicated in Table ||. 





Mass rang 


e (MeV/c 2 ) 


Nb. of 


Nb. of 


D sample 


M(D°/+) 


AM 


Df 


D£ + 


L> U -> K~TT + 


1820-1910 


>160 


752±41 


6±18 




1840-1890 


>160 


689±43 


39±26 


D + -> K-7T+1T+ 


1830-1910 


>160 


763±44 


66±19 


D*+ -> (K-ir + )iv+ 


1790-1940 


143.5-147.5 


416±24 


18±5 


D* + -> {K-u + iT + iT-)ixt 


1840-1890 


143.5-147.5 


303±21 


5±5 


D *+ (lf-7r+(7r ))7r+ 


1500-1700 


<155 


522±33 


15±12 



Table 1: Mass selections and number of D candidates observed in each decay channel 
(with their statistical error). Note that most of the D* + — > -D°7r+ candidates were 
removed from the selected D° sample; the D + sample also includes D* + — > D + tt° (or 7) 
decays. 



5 Semileptonic b decay rate into Dn£ X or D*n£ X 

In this section, a search for any Dtt^ final state is described, based on the impact 
parameter distribution of the n** candidates relative to the primary interaction vertex 
and using the D decay channels selected in the previous section. 

5.1 D**£~ selection 

The selection criteria for the additional 7T#* candidate were identical for all decay chan- 
nels. All charged particles with a momentum greater than 0.5 GeV/c and produced in the 
hemisphere defined by the D(ir^)£~ momentum vector were considered as 7r** candidates. 
The invariant D^^ir^i mass had to be smaller than 5.5 GeV/c 2 . The it** track was 
required to have at least 2 hits in the Vertex Detector. Its combined RICH and dE/dx 
information had not to be compatible with the kaon hypothesis. The impact parameter 
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of this 71"** relative to the previously fitted D(jr*)£ vertex was required to be smaller than 
100 /im. 

For each tt** candidate, the impact parameter relative to the primary interaction vertex 
was computed in the plane transverse to the beam axis. The sign of this impact parameter 
was defined with respect to the D{n*)£ direction. It was positive if the intercept between 
the tt** and the D(tt*)£ momentum vectors was downstream of the primary vertex along 
the D{ir*)£ direction, and negative if it was upstream JTT| . 

The tt** impact parameter distribution of simulated B semileptonic decays is shown in 
Figure 0a. Compared with charged particles produced in b quark fragmentation or gluon 
radiation in jets (see Figure ^o), tt** from B — > D**£~vi decays present a tail at large 
positive impact parameters due to the long B lifetime. 



5.2 Backgrounds 

For real data D**£ candidates, two sources of background had to be subtracted: 

• Fake D associated to a lepton candidate: this combinatorial background was es- 
timated by using events in the tails of the mass distributions of Figure [TJ, after a 
normalisation to the fraction of events below the D signals. Figures |]a-c present the 
impact parameter distributions of all pion candidates associated to a fake D (points 
with error bars) and to a D in the tails of the mass distributions (histograms) from 
the D samples selected in the qq simulation. A good agreement is found between 
the true background and the mass tail estimate. 

• True D associated to a fake lepton: this background is due to charged pions and 
kaons misidentified as leptons. It has been subtracted by using the it** candidates 
produced in the same direction as a wrong sign D£ + event (shown in the hatched 
histograms of Figure |l|) where the D candidate was selected in the mass range 
defined in Table [I]. Figure 0d presents the impact parameter distributions of all pion 
candidates associated to a fake lepton (points with error bars) and to a D£ + event 
(histograms) from the D samples selected in the qq simulation. Here also a good 
agreement is found between the true background and the £ + estimate. 

In the real data events, the same procedure was applied. The shapes of these backgrounds 
were taken from the real data themselves and their normalisation was estimated according 
to the fit of the mass distributions of Figure [I]. 

After the subtraction of these backgrounds, all the remaining pions can be attributed 
to b decays into Dtt£~X final state. However, four kinds of pions are still to be considered: 

• genuine tt** from B — > D**£~ue decays (see Figure ^|a); 

• particles from jet fragmentation (see Figure 0b); 

• "t, c — > £" background: it includes pions produced in D** decays when the D** is 
not issued from a direct semileptonic (e or //) b decay. This D** can be produced 
in b — > D**t~v t decay, or in b — > D**D^X (or b — > DD**X) transitions, when the 
other Dig-) (or D) meson decays semileptonically; 

• "hadronic" background: it is due to other hadrons, denoted H, produced from the c 
in b — > c — > £ decay events or from the c in b — > c — > £ (when the other charm quark 
fragments into a D meson). Such hadrons can be also emitted directly from the B 
meson. 

Despite the momentum and transverse momentum cuts applied to the lepton, these last 
two classes were not fully eliminated. Their impact parameter distributions were similar 
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to the impact parameter distribution of genuine tt^ from b semileptonic decays. These 
two last backgrounds were thus fitted together with the genuine tt** signal and subtracted 
only afterwards. Measured results were used for their rates and their selection efficiencies 
were obtained from the simulation (see Section 5.4.3). 

5.3 Total yield 

In the real data, the impact parameter distributions of the tt^ candidates of the 
"right" sign D°7t + £~, D + n~£~ and D* + n~£~ samples are shown in Figure |j. They were 
fitted, fixing the fake D and fake lepton backgrounds, but letting free the normalisation 
of the fragmentation and components. Figure || shows the same distributions, after 
subtraction of the fake D and fake lepton backgrounds. Similar fits were performed to 
the "wrong" sign D°tt~£~ , D + tt + £~ and D* + it + £~ samples and are shown in Figure |^ 
after the subtraction of the fake D and fake lepton backgrounds. 

Instead of rejecting kaons in order to select 7T#* , kaons were also identified in order 
to select from D sl -> D*°K + , D* 2 -> D°K + decays or any other D°K + X final state 
from other D s j resonances. The corresponding impact parameter distributions are shown 
in the same figures as above. The kaon rejection requirement led to a 7r** identification 
efficiency of (92 ± 1)% and a probability of wrong assignment as a kaon of (8 ± 1)%. The 
kaon identification requirement lead to a identification efficiency of (60 ±2)% and 
a probability of wrong assignment as a pion of (40 ±2)%. These factors were obtained 
from the real data, as explained in Section 5.4.1 (f n id correction). The numbers of fitted 
D°n£- and D°K£~ were 163 ± U(stat) and 39 ± 15(stat) (48 ± 21(stat) and 5 ± 8(stat)) 
in the "right" sign ("wrong" sign) samples. These values need to be corrected in order 
to take into account the fraction of kaons misidentified as pions. 

The final amounts, N(Dir£~), of "right" and "wrong" signs fitted candidates are pre- 
sented in Table ^| for all considered channels. In the D° channel, the separated contri- 
butions of D°tt, D°K and the total D°h are also indicated (where u h" means that the 
7r« candidate was selected without identification). Significant numbers of "right" sign 
candidates are fitted for all channels, whereas the number of "wrong" signs are clearly 
smaller. 



5.4 D**£ efficiency 

The semileptonic branching fraction of a b quark into Dn final state was then measured 
as follows: 

BR(6 - D«) = f f ±- f I^pi - T H - To, (2) 

1\Z £tib ^ ^Dl £** Jcor -DitD 

where Nz and ez are defined in Section 3, Rb = 0.2166 ± 0.0007 is the Z hadronic decay 
rate into bb events |T3]| ; the branching fractions, BRp, in the three decay modes BR(D° — > 
K-7c + ) = 0.0385±0.0009, BR(D + -> R-7i + 7i + ) = 0.090±0.006 and BR{D* + D°tt+) = 
0.683 ± 0.014 are used 0]. The efficiencies to reconstruct and select the D£~ and 7r« (or 
K^) candidates from B — > D**£~h>n decays, denoted ene and e** respectively, are indicated 
in Table ^| They were obtained from the simulation and corrected by the factors f cor 
and fr B which are described below. The correction factors f T . c ^e and Tu account for the 
"r, c — > £" and "hadronic" backgrounds introduced in Section 5.2; To* is the background 
due to residual D* + 7r~£~ which applies to the "wrong" sign D°n~£~ and D°K~£~ samples 
only. 
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Samnlp 




D u ti£- 


D U K£~ 


D + tt£~ 


D* + ir£~ 


1\ [ ngnt sign ) 


202 ± 37 


182 ± 39 


20 ± 18 


( ± ZD 


1 QO -I- OO 
L6Z ± ZZ 


J V ^ wi ong ^lgn I 


53 ±23 


55 ±24 


-2 ± 10 


A\ A- on 


OA _|_ i a 


£m 


0.127 ±0.002 


0.095 ± 0.002 


0.150 ± 0.002 


£** 


0.655 ±0.006 


0.649 ± 0.008 


0.654 ± 0.005 






1.87 ±0.09 




1 


3.02 ± 0.16 


P 

JM D 




0.94 ±0.01 




0.98 ± 0.01 


1.01 ± 0.01 


r 

JVD 




1.00 ±0.03 




i nn _i_ n no 
l.UU ± U.Uo 


1 nn _i_ n no 
l.UU ± U.Uo 


P 

JDvtx 




1 




r\ c\>~j i /~i no 

0.97 ± 0.03 


1 


f 

JKid 




0.84 ±0.02 




n sq -U n no 
U.oo ± U.UZ 


1 


firid 




1 




0.92 ±0.01 


U.92±U.01 


fr B ( "right" sign) 




1.02 ±0.02 




0.98 ±0.02 


0.98 ±0.01 


f TB ( "wrong" sign) 




0.98 ±0.02 




1.02 ±0.02 


1.03 ±0.03 




0.075 ±0.030 


0.075 ±0.030 


0.075 ±0.030 


Th (xlO" 3 ) 


1.06 ±0.29 


0.78 ±0.22 


0.28 ±0.09 


0.40 ±0.12 


0.41 ±0.11 


^ (xlO" 3 ) 


0.25 ±0.06 


0.23 ±0.06 


0.02 ±0.01 









Table 2: Number of fitted Dtt£~ candidates; reconstruction times selection efficiencies of 
the D£~ and tt** (or K**) from B —> D**£~vi decays; correction factors introduced in 
equation (2). Errors are statistical only (except for f T>c -^i and Th)- Note that most of the 
D* + — > D°tt^ candidates have been removed from the selected D° sample; the D° sample 
also includes D*° —> D°tt (or 7) decays; the D + sample also includes D* + — > D + tt° (or 
7) decays. 

5.4.1 Efficiency correction 

The correction to the reconstruction and selection efficiency is expressed as 

/cor = fx-ir fM D fvD fDvtx /if id firid'- 

• In the D°£ and D* + £ samples, only the Kir decay channel was used to estimate the 
em efficiency For these samples, fxn = N(D£)/N(D Kw £) where 

N(D£) = N(D£-) - N(D£ + ) is the difference between the total number of D£ 
candidates quoted in Table [I] and N(D Kn £) is the same difference computed in the 
Kit decay mode only of the D°. In the D + £ sample, fxir = 1- 

• Due to the D mass ranges required in Table [I], fu u accounts for the mass width 
differences observed in real data and simulation. 

• A large sample of D* + — > (K~ti + )'k^ reconstructed inside 6-tagged jets was used in 
order to estimate selection efficiencies related to the detector response: the Vertex 
Detector information which was required for all channels (/vd factor), the vertex 
quality cuts for the D + sample (fDvtx), and the kaon identification for the D° and 
D + samples (fxid)- For the study of the K^n + n + vertex quality in the D + sample, a 
three tracks K^-k + ti^ vertex was also fitted in the dedicated D* + sample and similar 
cuts were applied. 

• The kaon rejection (or identification) requirement of the 7r** (-K**) candidates was 
also checked on the same dedicated D* + sample and a correction factor f^ was 
inferred. 
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5.4.2 B lifetime correction 

The difference between the known values of the B mesons' mean lifetimes (t b + = 
1.65 ± 0.04 ps, r B o = 1.56 ± 0.04 ps |7|, t b Jt b o ~ 0.99 - 1.01 and that used in the 
Monte Carlo simulation (V|f c = 1.6 ps) has two consequences: 

• It affects the decay length selection efficiencies described in Section 4. But as these 
selections were applied to the decay lengths divided by their errors, the relative 
correction to t B n was found to be of about ±0.2% only. It was thus included in the 
following f Tg factor. 

• It also affects the shape of the impact parameter distribution of simulated 
n„ candidates which is used to fit the amount of Dtt£~ candidates in real data. 
The distribution shown in Figure 0a was thus recomputed by weighting each simu- 
lated event and by using the B° mean lifetime for "right" sign D°h + £~ and "wrong" 
sign D^ + ir + £~ candidates, or the B + mean lifetime for "right" sign D^ + n~£~ and 
"wrong" sign D°h~£~ . These new 7T*# impact parameter shapes were used to fit the 
real data distributions shown in Figures |4]-|5|. The difference between the number of 
fitted D^7f„£X candidates observed with and without the weighting procedure is 
described by the correction factor, f TB , given in Table |2|. 



5.4.3 Physical background correction 

The physical background contributions are determined in the following way: 

• According to the simulation, still (7.6 ± 0.4 (stat))% of D* + — > -D°7r+ remained in 
the D° sample; this value was used to determine the JF D , factor. 

• The fraction of b — > r — > £ events is evaluated as: 

fr^t = " " ' v , — — = 0.0075 ± 0.0020 3 

BR(o — > £ h>eX) €£ 

where BR(b -> t~v t X) = (2.6 ± 0.4)%, BR(r" -> £~v t v T ) = (17.64 ± 0.06)% and 
BR(6 -> t~v t X) = (10.99 ± 0.23)% 0; e T ^ e /e e = 0.18 ± 0.04 (stat) is the ratio of 
the lepton selection efficiencies in b — > t~v t X and b — > l~vgX simulated events. 
The fraction of b — > c — > £ events is evaluated as: 

BR(6^c^£) e c ^ £ 

where BR{b -> c -» £)) = (1.6 ± 0.4)% @; e^/e/ = 0.32 ± 0.02 (stat) is the ratio 
of the lepton selection efficiencies in b —> c — > £ and b — > £~V(X simulated events. 
The fraction of b — > c — > ? events is obtained as previously, but the probability, 
Pvk-»z), for the virtual W 7- to decay into a D oy D~ meson has to be taken into 
account: 

f~t = m%~*f~n Pw - D — = °-° 20 ± °-°° 5 (5) 
BR(o — > £ i^X) q 



where BR(6 -> c ->• £) = (7.8 ± 0.6)% and = (9.0 ± 1.9)% H. Finally 



the fraction of P** not issued from a direct semileptonic b decay is evaluated to be: 
f T ^ t = f T ^ e + h-*l + fc->e = 0.075 ± 0.017 ± 0.025 . (6) 

In the first error, which is the sum of the uncertainties quoted in equations (3-5), 
the errors on fc^e and / c _>^ have been added linearly because Pw-*d was used in 



10 



reference [14 to evaluate BR(5 — > c — > £). The second error in equation (6) is an 
estimate of the uncertainty due to possible phase space or QCD corrections between 
the b — > £ and the b — > r, c — > £ decay channels with a D** in the final state. 
The "hadronic" background is evaluated as: 

T H = BR(6 -> DX) (1 - P^ D ) BR(6 -> (c or c) -> £±) e ^ Dmx (7) 



with 



BR(6 -> (c or c) -> £±) = BR(6 -> c -» £) + BR(6 -> c -> £) P, 



W->D 



(2.30 ±0.56)% 



BR(6 -» P°X) = (60.1 ± 3.2)% Q, BR(6 -> P+X) = (23.0 ± 2.1)% and 
BR(6 -> P*+X) = (23.1 ± 1.3)% pi; the difference BR(6 -> P°X)-BR(6 -> 
P*+X)BR(P*+ -> P°7T + ) = (44.3 ± 3.3)% is used for the P°r analysis, where the 
P* + were rejected. In the Dti£~ analyses, €b~,DHex/(^De^**) = 0.084 ± 0.010 (stat) 
was determined from the simulation as the ratio of selection efficiencies be- 
tween hadrons from charm decay in b — > c — > £ events, and genuine in 
P — >• D**£~U£ decay; in the D°K£~ analysis, this ratio was estimated to be 
0.030 ± 0.006 (stat). The resulting Tn values are reported in Table 0. 



5.5 Systematics 

The systematics are detailed in Table |. As a cross-check of the procedure, the same 
analysis was repeated on simulated qq and bb samples: 

• 1998 ± 107 (1017 ± 73, 870 ± 62) D** candidates were fitted in the "right" sign 
D°7i + £- (D + 7i-£-, D*+tt-£~) samples whereas 1934 (1106, 879) D** were expected; 

• 396 ±67 (219 ±52, 60 ±38) P** candidates were fitted in the "wrong" sign D°ti + £- 
(D + ti-£-, D*+n-£-) samples whereas 333 (235, 62) P** were expected. 

A good agreement was thus obtained in the simulation between the fitted and expected 
contributions, the related statistical error being used to estimate the systematic 
uncertainty due to the subtraction of the fake P and fake lepton backgrounds. The 
remaining statistical error of the Monte Carlo simulation is due to the limited number of 
generated P — > Diri~p^X events. 

Following the detailed study of reference |TB|, a ±0.3% uncertainty is assigned to the 



reconstruction efficiency of each charged particle. 

The uncertainty on the impact parameter resolution has two sources: 

• Impact parameter relative to the primary interaction vertex: 

the uncertainty on the parameters a and b of equation (1) affects the impact pa- 
rameter distributions of Figure and thus the result of the fit to the real data; the 
corresponding relative systematic error is estimated to be at most of ±1%. 

• Selection of the impact parameter of the 7r** candidate relative to the D(7r*)£~ vertex: 
this impact parameter was required to be smaller than 100 /im, which allowed the 
selection of about 82% of genuine 7r** candidates. A variation of ±10% of the 
impact parameter resolution modified the relative efficiency by about ±2.5%, the 
magnitude depending on the considered P decay channel. A similar uncertainty 
was also inferred by comparing, in D* + £~ real data and simulation, the effect of a 
cut on the impact parameter of the lepton relative to the D°7r*£~ vertex. 
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The overall systematic uncertainty due to the impact parameter resolution is given in 
Table |3|. 

The uncertainty due to the 7r** momentum spectrum was evaluated by comparing the 
7r« selection efficiencies in simulated Dj — > Dtc and Dtttt decays: a relative difference of 
±1% was observed. 



Error source 


D U M~ 


D l) Tx£r 


D U K£- 


D+7i£- 


D* + 7c£~ 


b — > £ decay model 


13| 


±1.2 


±1.2 


±1.2 


±1.2 


±1.2 


TB |17| 




±1.9 


±1.8 


±2.2 


±1.6 


±1.2 


7T** momentum spectrum 


±1.0 


±1.0 


±1.0 


±1.0 


±1.0 


r, c — > £ background 


±3.5 


±3.5 


±4.2 


±3.5 


±3.5 


"hadronic" background 


±2.5 


±2.1 


±9.2 


±2.4 


±2.3 


BR(D° -> 


i 


±2.3 


±2.3 


±2.3 


— 


±2.3 


BR(D + -> K-U+U+ 


) E 




- 




±6.7 


- 


BR(D* + D°7T + ) 


E 










±2.0 


fake Z)£ backgrounds 


±5.4 


±5.4 


±5.4 


±7.2 


±7.1 


MC statistics 




±2.0 


±2.0 


±2.3 


±2.6 


±1.7 


track reconstruction 


±1.3 


±1.3 


±1.5 


±1.6 


±1.6 


impact parameter resolution 


±3.2 


±3.1 


±3.7 


±3.7 


±2.8 


mass resolution 




±1.2 


±1.1 


±1.4 


±1.1 


±1.1 


VD requirement 




±3.3 


±3.2 


±3.9 


±3.2 


±3.3 


D vertex selection 








±3.3 




K (from D) identification 


±2.6 


±2.6 


±3.1 


±2.6 




7r« (or K**) identification 




±2.1 


±18.1 


±1.2 


±1.2 


lepton identification 


±2.0 


±1.9 


±2.3 


±2.0 


±2.0 


Total 


±9.9 


±9.9 


±22.9 


±13.3 


±10.6 



Table 3: Relative systematic uncertainties (%) on the b semileptonic branching fractions 
into D**£~ui final states ("right" sign only). 



5.6 Results 

From the previous study, the b semileptonic branching fraction can be computed in 
each Dir£~ or D°K£~ final state. The corresponding results are reported in Table [| which 
includes the statistical and systematic errors. The "right" sign values are in agreement 
with those measured by the ALEPH collaboration ||, except for the D°tt + £~X channel 
where the DELPHI result is two standard deviations larger. 

The "wrong" sign results are at less than 2 standard deviations from zero, thus Ditn 
final states will be neglected in the following. The D°K£~ production rate is also found 
to be compatible with zero. Thus only Dn£~ final states will be considered in the follow- 
ing. As a further cross-check, Tables [5|-[| present the b semileptonic branching fraction 
measurement for electrons and muons separately and for the various D decay channels. 

Using the production fraction BR(b -» B°) =BR(b -> B~) = 0.395 ± 0.014 0, the 
following branching fractions are measured: 

BR(B° DVr^) + BR(B° D*°ir + £-i7 e ) = (2.70 ± 0.64 (stat) ± 0.28 (syst))% 
BR{B- -> D + n-£~i7 e ) + BR(S" -> ZT+Tr^r^) = (2.08 ± 0.47 (stat) ± 0.20 (syst))% 
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BR(b -> D**£~V 
D°h£~X 


{ ) (xl0~ a ) 
D°ni-X 


D°K£~X 


D+n£~X 


D*+ir£-X 


DELPHI "right' 
11.6 ± 2.4 ± 1.1 


sign 

10.7 ±2.5 ±1.1 


1.0 ± 1.1 ±0.2 


4.9 ± 1.8 ±0.7|4.8 ±0.9 ±0.5 


DELPHI "wrorif 
1.9 ± 1.4 ±0.4 


f sign 
2.3 ±1.5 ±0.4 


-0.4 ±0.6 ±0.1 


2.6 ± 1.5 ±0.4 


0.6 ±0.7 ±0.2 


ALEPH "right" 


sign 
4.7 ± 1.3 ± 1.0 


2.6 ±1.2 ±0.8 


3.0 ±0.7 ±0.5 


4.7 ±0.8 ±0.6 



Table 4: Semileptonic branching fractions BR(b — > D**£~i/z) measured in DELPHI for 
each Dir£~ or D°K£~ final state. Similar results obtained in ALEPH are also pre- 
sented [0. The first errors are statistical and the second systematic. Note that D°s 
from the D* + — > D°tt+ decay mode are removed from the D°£~ results, which still in- 
clude D°s from D*° — > D°ti (or 7) decays; results on D + also include D* + — > D + ti (or 
7) decays. 



BR(b -> D**£~P £ ) (xl0~ a ) 


D"h + £-X 


D + n-£-X 


D* + n-£-X 


e 


8.8 ±3.0 


5.3 ±2.3 


5.1 ± 1.2 




14.6 ±3.4 


4.5 ±2.5 


4.4 ±1.2 


Average 


11.6 ±2.4 


4.9 ± 1.8 


4.8 ±0.9 



Table 5: Semileptonic branching fraction for electrons and muons separately. Errors are 
statistical only. 



BR(b ^ D**£-p e ) (xl0~ a ) 


D u h + £~X 


D* + 7i-£-X 


Kn 

K7T7T7T 

Knn 


13.6 ±2.6 
9.8 ±4.1 


5.5 ±1.5 
5.7 ± 1.8 
3.9 ±1.3 


Average 


11.6 ±2.4 


4.8 ±0.9 



Table 6: Semileptonic branching fraction for the different D decay channels. Errors are 
statistical only. 



According to isospin conservation rules and assuming that only Dn and D*ir final 
states contribute, the ratios between final states involving charged and neutral pions are 
predicted to be: 



D°7l + + D*°7T + D + 7T- + D* + 7l- 



(9) 



L> 7T + D 71 D i\ + D T\ 
allowing the following branching fractions to be inferred: 

BR(B° -> Dti£-u £ ) + BR(£?° -> D m ir£~j7 t ) = (4.05 ± 0.96 (stat) ± 0.42 (syst))% (10) 
BR(B~ -> Dn£-i7 e ) + BR(B~ -»- DWi/i) = (3.12 ± 0.71 (stat) ± 0.30 {syst))% (11) 
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These values are in good agreement. Neglecting Dirir final states and using as a constraint 
the equality of equations (10) and (11), the overall B meson semileptonic branching 
fraction into any D^ir final state can be obtained: 

BR(B D-kI-vi) + BR{B -> D*ir£~i7 t ) = (3.40 ± 0.52 (stat) ± 0.32 (syst))% . 

Assuming that the isospin invariance used in equation (9) applies also to Dn and D*n 
separately, the following branching fractions are also inferred: 

BR(B -»■ Dirl'vi) = (1.54 ±0.61 (stat + syst))% 
BR(B -»> D*Tt£~u e ) = (1.86 ±0.38 (stat + syst))% 

with a correlation coefficient of -0.33 between the results. 



6 Overall b decay rate into D£ v^X final states 

In this section, a measurement of the semileptonic branching fractions BR{b — > D£~X), 
where D stands for D°, D + or D* + , is presented. The method used is similar to that 
described in Section 5.4: 

BR(6 DtX) = f - J- -J- Izip£ _ ^, + ^ (12) 

where N(Dk-k£) = N(D£~) — N(D£ + ) is the difference between the total number of D£ 
candidates quoted in Table [I], using only the Kit (Ktttt, Kittt*) decay mode in the D° 
(D + , D* + ) analyses; f' cor = f Mo f VD f Dvtx f Kid ; T' D * is only used for the D° sample where 
a fraction of (7.6 ± 0.4 (stat))% of D* + — >• D°n+ decays is included. 

According to the simulation, the reconstruction efficiency of D£ final states depends 
slightly on whether or not the D meson originates from a D**. In the absence of D**, 
the reconstruction times selection efficiency, eoe, of Table |2] has to be multiplied by the 
factor r] De = 1.08 ± 0.02 (1.13 ± 0.02, 1.07 ± 0.02) for a D° (D+, D* + ) final state. Thus 
the observed production fraction of Dtt£~ and D°K£~ final states (denoted Dh£~) has 
to be taken into account in equation (12) and the following factor is introduced: 

J~dm = T]m ^ 1 BR(b -> DhtTX) . (13) 

VD£ 

The overall b semileptonic branching fractions are thus measured to be: 

BR(fe -> D°£-p t X) = (7.04 ± 0.34 (stat) ± 0.36 (syst.exp) ± 0.17 (BR D ))% 
BR(b -> D + e~9 e X) = (2.72 ± 0.19 {stat) ± 0.16 (syst.exp) ± 0.18 (BR D ))% (14) 
BR(b -> D* + £-9iX) = (2.75 ± 0.17 (stat) ± 0.13 (syst.exp) ± 0.09 (BR D ))% 

where the D°£~ result includes also D° coming from D*° and also (contrarily to Section 5) 
D* + — > D°7r^ decays, the D + £~ result includes also D + coming from D* + — ► D + tt° (or 7) 
decays and A means "anything" (possibly a hadron coming from a D**). These results 
are compared in Table [7| with those measured by the OPAL collaboration [0: the D° and 



D* + values are in agreement whereas the D + results present a difference of two standard 
deviations. The systematics are detailed in Table [8|. 

The relative yield of D* + i{£~b>iX over all D* + £~i?fX is a contribution to the systematic 
uncertainty of various measurements, particularly of V c b PflpOl. From Table [| and 
equations (9) and (14), the following ratio is obtained: 

BR(6 -> D* + tt-£-u £ X) + BR(b D* + 7r°r^A) 



BR(6 -> D* + rv t X) 



0.26 ± 0.05 (stat) ± 0.02 (syst) 
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which significantly improves on a previous DELPHI measurement 0. 



BR(b -> D£'P e X) 
(%) 


DH~ 


D+£~ 


D*+£~ 


DELPHI 
OPAL 


7.04 ±0.34 ±0.36 ±0.17 
6.55 ±0.36 ±0.44 ±0.15 


2.72 ±0.19 ±0.16 ±0.18 
2.02 ±0.22 ±0.13 ±0.14 


2.75 ±0.17 ±0.13 ±0.09 
2.86 ±0.18 ±0.21 ±0.09 



Table 7: Overall semileptonic branching fractions into Dl~ final states as measured in 
DELPHI and OPAL fl9fl . The first errors are statistical, the second are experimental 
systematics and the last are due to the exclusive D branching fractions, BR D . 



Error source 


dH- 


D+£~ 


D*+£- 


b — > £ decay model 
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±1.2 


±1.2 


±1.2 


TB 




±0.2 


±0.2 


±0.2 


r, c — > £ background 


±1.8 


±1.8 


±1.8 


BR(D° -> K-n + ) 


i 


±2.3 




±2.3 


BR(D + -» A'-7T + 7T + 


) 1 




±6.7 




BR(L>* + -> D°7T + ) 


I 






±2.0 


MC statistics 




±2.4 


±2.4 


±2.7 


track reconstruction 


±0.9 


±1.2 


±1.2 


mass resolution 




±1.1 


±1.0 


±1.0 


VD requirement 




±2.0 


±2.0 


±2.0 


D vertex selection 




±3.1 




K identification 




±2.4 


±2.4 




lepton identification 


±1.8 


±1.8 


±1.8 


Total 


±5.5 


±9.0 


±5.6 



Table 8: Relative systematic uncertainties (%) on the b semileptonic branching fractions 
into D£~i?iX final states. 



7 Summary and conclusion 

Using DELPHI data recorded from 1992 to 1995, the overall b semileptonic branching 
fractions into D°, D + or D* + final states have been obtained: 

BR(b -> D°rviX) = (7.04 ± 0.34 (stat) ± 0.36 (syst.exp) ± 0.17 (BR D ))% 
BR(6 -> D + £-9tX) = (2.72 ± 0.19 (stat) ± 0.16 (syst.exp) ± 0.18 (BR D ))% 
BR(6 -> D* + £-9tX) = (2.75 ± 0.17 (stat) ± 0.13 (syst.exp) ± 0.09 (BR D ))% 

where the D° and D + results include also contributions from D*° and D* + decays. 

Evaluating the yield of charged pions from higher excited states or from non-resonant 
D^n final states, the following branching fractions have been measured: 

BR(6 -> D°n + £-i7 e X) = (10.7 ± 2.5 (stat) ± 1.1 (syst))l()- 3 
BR(b -> D + ir-£-p e X) = (4.9 ± 1.8 (stat) ± 0.7 (syst))l()- 3 
BR(b -> D* + 7r^rz7 £ X) = (4.8 ± 0.9 (stat) ± 0.5 (syst))10~ 3 
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and 

BR(b -> D^Tt-TutX) = (2.3 ± 1.5 (stat) ± 0.4 (s?/st))10~ 3 

BR(6 -> D + ir + ri7 e X) = (2.6 ± 1.5 (stat) ± 0.4 (syst))10" 3 

BR(6 D* + 7T + £~i7 e X) = (0.6 ± 0.7 (stat) ± 0.2 (s?/st))10~ 3 

where the _D* + — > D°ir^ decay mode is not included in the BR(fe — > D°-k l~i>iX) results. 
Neglecting D7T7T final states and assuming isospin invariance, the separated branching 
fractions are inferred: 

BR(B -> DTxt~v t ) = (1.54 ± 0.61 (stat + s?/st))% 
BR(B -> rritfvt) = (1.86 ±0.38 (stat + syst))% . 

The measured overall branching fraction: 

BR(B -> Diri've) + BR(B -> D*nri7 e ) = (3.40 ± 0.52 (stat) ± 0.32 (syst))% 

is found, in good agreement with the expectation from the difference [|TJ: 

BR(5 -> r£pT) - BR(5° -> D+ru e ) - BR(B° -> D* + rz7 £ ) = (3.85 ± 0.42)% 

but is larger than a previous ALEPH result of BR(B — > D7r£ _ z/ £ ) +BR(_B — > D*n£~i7 e ) = 
(2.26 ± 0.29 (stat) ± 0.33 (syst))% @. 
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figure 1: Invariant mass distributions in the (a) D° — > K 7r + (b) D° — > A' 7r + 7r + 7r 
(c) L> + — > A' _ 7r + 7r + and (e) _D* + — > (fr~7r + 7r + 7r~)7r+ decay channels; mass differ- 
ence distributions M(K-tt + tt+) - M(K-n + ) in the (d) D* + -> (7^+)^+ and (f) 
D* + — > (A' _ 7T + (7r ))7r+ decay channels. The reconstructed .D* + candidates have been 
removed in a,b,c. Right charge D£~ (dots) and wrong charge D£ + (hatched histogram) 
events are shown. The solid line curves are fits which include a background parameteri- 
sation (dashed curve alone) and Gaussian functions for the signal (see Section 4). 
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a) MC 
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figure 2: Impact parameter relative to the primary interaction vertex in simulated B 
semileptonic decays for a) 7r« from D** decay (using a B mean lifetime value of 1.6 ps) 
and b) charged particles from jet fragmentation (see Section 5.1). 
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figure 3: In the simulation: impact parameter distributions of pions accompanying (a-c) a 
fake D meson (points with error bars) oraD selected in the tails of the mass distributions 
(histograms); (d) a fake lepton (points with error bars) or a D£ + where the D was selected 
in the signal region (histograms, see Section 5.2). 
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figure 4: Impact parameter relative to the primary interaction vertex in real data for 
"right" sign D°n + , D°K + , D + n~ and D* + n~ candidates. The black and cross-hatched 
histograms are the estimated contributions from fake leptons and fake D mesons, re- 
spectively. The hatched and empty area histograms are the fitted contributions from jet 
fragmentation and 7r« from D** decays, respectively (see Section 5.3). 
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figure 5: Impact parameter relative to the primary interaction vertex in real data for 
background subtracted "right" sign D°n + , D°K + , D + n~ and D* + n~ candidates. The 
hatched and empty area histograms are the fitted contributions from jet fragmentation 
and from D** decays, respectively (see Section 5.3). 
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figure 6: Same as Figure [| for background subtracted "wrong" sign D°ti , D°K~ 
D + n + and D* + n + candidates. 



